Competing order parameters and a tricritical point 
with a difference 



Arghya Dutta*, J. K. Bhattacharjee 

S. N. Bose National Centre for Basic Sciences 
Block- JD, Sector-3, Salt Lake, Kolkata, Pin-700098, India. 



Abstract 

We propose a mean-field, phenomenological Ginzburg-Landau free energy func- 
tional with two competing order parameters for a two-component, spin-polarized 
Fermi gas. This free energy supports a tricritical point which is different from 
the conventional one and this change offers a correct understanding of the ex- 
perimental phase diagram of imbalanced Fermi systems ( Shin et al, Nature, 
2008). The specific heat also happens to be different than in standard theory. 
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1. Introduction 

Population imbalance between spin-up and spin-down fermions in a mixture 
of two Fermi gases leads to exotic superfluid phases. In past ten years, these 
exotic superfluid phases have been extensively studied, both theoretically 0, 0, 
S 0, HI and experimentally 0, 0, @, S E3, El in ultra-cold atom ic g ases and in 
quark-gluon plasma which resides in the core of neutron stars[l2l Il3j|. Simi- 
lar imbalanced fermionic systems were studied in electron superconductors in a 



magnetic field as early as nineteen-sixties[14j, [15|, and showed the possibility of 
tricritical points where the second-order and the first-order lines meet along with 
the line separating stable and unstable superconducting phases. That the tri- 
critical point is fundamental to the understanding of superfluidity of polarized, 
two-component Fermi gases was first pointed out by Parish et al[16||. 

The standard model for studying tricritical point starts with the following 
Ginzburg-Landau free energy per unit volume: F — ^ip 2 + jtp 4 + §?/> 6 , where 
ip is superfluid order parameter, a is scaled temperature, c is a positive con- 
stant, and b can take both positive and negative values. Near the tricritical 



'Corresponding author. Tel: +91 (033) 2335 5706 fax: +91 (033) 2335 3477 
Email address: arghyaObose . res . in (Arghya Dutta ) 



Preprint submitted to Elsevier 



July 17, 2012 



Unstable 
(Coloured 
rag ion) 



b 




Normal 



11 



a 



Stable 
Superfluid 



Figure 1: Standard 0(»/> 6 ) theory phase diagram as a function of the scaled temperature (a) 
and temperature-independent, polarization-dependent fourth order coupling constant(b). The 
shaded region represents a thermodynamically-unstable, finite-i/j region. 

temperature (To), one usually expands a as a power series of temperature dif- 
ference - a = ciq(T — To) where ao is a constant and T is the temperature, 
and then searches for extrema in the free energy landscape. If all the coeffi- 
cients of the free energy become positive then F attains its minimum value at 
ip = 0. In other words, the system remains in the normal state. For b < 0, 
however, the extrema of the free energy are located at those ^-values which are 
solutions to the equation a + bip 2 + ctp 4 = 0, specifically the non-zero minima 
are located at tp 2 = —b/2c ± \/b 2 — 4ac/2c. Clearly a non-zero minimum ex- 
ists if b 2 > Aac. The minimum will be unstable energetically unless it satisfies 
= F(ip = 0) = F(ip 2 = -b/2c+Vb 2 - 4ac/2c), which happens if ac = 36 2 /16. 
So, in summary, in the free energy parameter space, we have an unstable finite- 
ip phase between a — 3b 2 /16c and a = b 2 /4c. The corresponding phase diagram 
is plotted in fig. dl]) 

Recently, an experiment by Shin et al[17] captured the first-order transition 
very clearly by measuring spatial discontinuity in spin polarization . In the 
experimental phase diagram reported by Shin et al [see fig. ([2])], we note that 
the width of the unstable, finite-^ region increases with decreasing temperature. 
It is, however, apparent that whatever polarization dependence we assign to the 
temperature-independent constant b of the standard 0(ip 6 ) theory, the width of 
the unstable superfluid region decreases with decreasing temperature. 

Keeping this discussion in mind, we set out to reconstruct the free energy 
of this imbalanced Fermi system following the work done on antagonistic or- 
der parameter in superconductor-ferromagnet phase transitions, several decades 
ago by Blount et aljlij]. Our phenomenological model free energy is a general 
Ginzburg-Landau type free energy with two order parameters, which are actu- 
ally embedded in the system: one superconducting(^), and another imbalance 
parameter (to). This free energy explains the phase diagram properly and offers 
some valuable physical insights about imbalanced Fermi systems. 
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Figure 2: Relevant portion of the phase diagram of spin-imbalanced, two-component Fermi 
gas, as reported by Shin et al(Nature,2008). The point 'T' denotes the tricritical point. The 
most important characteristics, which we are trying to point out in this paper, is that the 
width of the unstable region (shaded violet) increases with decreasing temperature, whereas 
the width of the unstable region, as obtained from a standard 0(ip e ) theory, decreases with 
decreasing temperature, as can be seen in fig.JTJ. 

2. The Model and Results 

Our model free energy per unit volume is 

F = -ib 2 + -ib 4 + — m 2 - mh + — m 2 ib 2 (I) 

in which if) is superfluid order parameter, a is given by a — a Q (T — Tq) where ao 
is a constant and Tq is the normal-supcrfluid transition temperature in absence 
of any imbalance. The parameter b is a temperature independent constant. 
As usual in mean field theory, we have a uniform superfluid state for T < Tq 
with ip 2 = —a/b. The imbalance parameter(m) is the difference in the number 
of spin-up and spin-down fcrmions. This parameter represents magnetization 
in the superconducting scenario and polarization in the ultra-cold atomic gas. 
Following the standard Ginzburg-Landau formalism, we include a term propor- 
tional to m 2 in the free energy (where the coefficient A is a constant) and the 
magnetic field (here chemical potential difference between the fermionic species) 
gives rise to an additional contribution —mh. We have neglected 0(m 4 ) terms 
in the free energy as high imbalance destroys superfluidity. Finally, the m 2 ip 2 
term, proposed after considering the symmetries of the system, signifies the in- 
teraction of the imbalance parameter with the superfluid order parameter. One 
has this freedom of choosing particular form of interaction while one is trying to 
explain the phase diagram of an imbalanced, two-species Fermi system at uni- 
tarity, because the particular functional form of the interaction does not matter 
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Figure 3: This figure depicts our model GL free energy(F) as a function of ip and m. This 
clearly shows the double minima of F with respect to tp and one minimum for m. 



in this resonantly-interacting regime. We have plotted the free energy in fig. 
©• 

To derive results for this model, we first minimise the Ginznurg-Landau 
free energy. Minimization results in two conditions on the order parameter - 
m = h/(A + Bip 2 ) and tp(a + bip 2 + Brh 2 ) = 0. Now for superfludity to be 
present in the system, one needs to solve the above equations simultaneously 
for non-zero solutions of the superfluid order parameter, and that is equivalent 
to solving 

(a + bip 2 )(A + Bip 2 ) 2 +Bh 2 = 0. (2) 

This result could also be viewed in a different way. We could integrate out the 
to- variable from Eq. (fl}, leading to an effective free energy F e ff- Minimising 
this with respect to ip again leads to Eq. This shows that results from 

conventional tricritical point can differ when ip 2 is not small and our contention 
is to show that this is the point of difference from standard tricritical point. 

From eq. ([2]), we note that the superfluid order parameter vanishes con- 
tinuously along the curve A 2 a + Bh 2 = 0. For temperatures greater than the 
tricritical temperature, a second order phase transition takes place as one crosses 
this curve to go from a normal(?/> = 0) to a superfluid phase^ ^ 0) or vice versa. 

As our primary aim is to explain the phase diagram found by Shin et al, we 
will now check whether there exists a first order normal-superfluid transition in 
this system. If there exists a non-zero solution of the superfluid order parameter 



See appendix for detail calculations. 
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Figure 4: Spin polarization(m) vs. temperature(i = T/Tq) phase diagram for a imbalanced 
two-component Fermi gas with resonant interactions. The area enclosed between the purple 
and blue line represents an energetically unstable region. The red dot is the tricritical point. 
For temperatures greater than the tricritical temperature superfluid-normal phase transition is 
second-order, while below it superfluid-normal phase transition is first-order which encounters 
an unstable phase en route. (Plotted using A=1.0, B=0.5 and b=0.2) 



in the equation 

F(i/>,m) =F(ip = a,m\^= ), (3) 

there will be a first-order phase transition in the system. We insert our particular 
free energy in this equation and found the following values for the superfluid 
order parameter 




Using the functional form of ^ 2 , we get two equations which relate the 
magnetic field and critical temperature : h 2 — ~aA 2 /B and h 2 = A(bA — 
2aB) 2 /8B 3 , across which a first-order normal-superfluid transition takes place. 
One of equations for the first order line is identical with the second-order line. 
The other curve is the purple line shown in fig®. The tricritical temperature (To) 
can be easily determined by evaluating the intersection point of the two branches, 



which works out to be ao = A 



[jb ~ l ) ± V 1 _ B 



We designate this point 



as a tricritical point because it is the end point of a series of critical points, i.e., 
the order of the phase transition changes from second to first at this point. 

Now let us figure out if there is any unstable region in the phase diagram 
within the scope of our model. This can be found out in this way- as we are 



5 



treating this model at a mean field level, the phase of the superfluid order 
parameter is unimportant. Hence by demanding the square of the superfluid 
order parameter to be real, we find that there exist a unstable state for a range of 
temperature and magnetic field. According to the experimental data of Shin et 
al, in this region phase separation takes place. But from our free energy, which 
does not include any gradient term, we can only conclude that this region is 
unstable. 

Another interesting feature of this model Ginzburg-Landau free energy is its 
specific heat: we get a rather peculiar behaviour of this quantity near tricritical 
point. Generally in 0(ip 6 ) theories, where tricritical point is quite generic, at 
temperatures just below tricritical point specific heat diverges as C ~a"5. 

To calculate the specific heat of this our system, we differentiate twice the 
effective free energy (F e ff) with respect to a, which is nothing but scaled tem- 
perature, and we get 

c= d 2 F eff = p] + «WUW 2 + ^/3V 



(5) 



da 2 \ da J 2 \ da 2 J 2 \ da J 2 \ da 2 

B 2 h 2 (d^ 2 \ 2 Bh 2 (d 2 i? 



(A + Bip 2 ) \ da J 2(A + BiP 2 ) 2 \ da 



So, as we can see from above equation, the nature of the specific heat is actually 
determined by the dependence of ip 2 on a which can be extracted from eq. (j4]). 
There are two distinct and interesting cases where the specific heat behaves 
differently. 

Near transition point of the superfluid phase diagram, ip 2 varies as a 1 / 2 , in 
general, which can be derived from eq. Q. If we put this a-dependence of ip 2 
in eq. (0, we get C=l+|+^— (A+Ba) 3 > wmcn : near the transition point, 
becomes 

So the specific heat have a jump discontinuity at the transition point. We em- 
phasise, this is different from the standard tricritical behaviour of specific heat, 
which actually diverges as we approach the tricritical point from lower temper- 
atures. 

Now let us consider the second case. Interestingly, we can expand our free 
energy to get a 0(tp 6 ) theory (for the steps see appendix). The coefficient of the 
ip 4 term will be equal zero, which is required for a standard tricritical point, if 
hi = |^r(from eq. (|A.3|) V Now ip will be proportional to a 1 / 4 , like in standard 
tricritical theory, if we set = -g + e, where e is a small number and in that 
case ip 2 reduces to 

12 e I / " 77 h 2 A 2 7I\ 



Putting the value of in the above equation, we obtain ip 2 = 4 + k \/£ 2 + 



G 



As e will tend to zero, if 2 wm diverge as y/e: finally we get the tp ~ a 1 / 4 de- 
pendence. Doing a standard tricritical analysis we getC~U+^- , as 

expected. So indeed there is a point is the phase diagram, whose (a, h) coordi- 



nates are I — =^r, y §§2 J , where we recover the standard tricritical behaviour. 

Thus specific heat measurement will thus serve as an interesting check on 
whether the tricritical point is of the conventional type(O(-0 6 ) theory) or the 
findings over here. The main point of the last few paragraphs is that our model 
free energy encompasses the standard tricritical phenomena and, at the same 
time, points to some interesting physics which standard tricritical free energy 
does not account for. Before proceeding to a comparison of our free energy with 
the available microscopic calculation of this systems, we want to make one final 
comment about the special case in which we set b — B in our free energy In 
this case, the transition point of our model and the standard tricritical point 
becomes one, and one can have C ~ a~ 5, as in standard tricritical case, also 
for the transition point of our model. 

To make the results quantitative, we connect the proposed free energy to 
the available microscopic calculations of the system [15j. The microscopic cal- 
culation was done for a superconductor with an internal magnetic field gen- 
erated by dilute magnetic impurity ions , neglecting the orbital effects, which 
was a justified assumption in a two-dimensional superconductor. The magnetic 
field differentiates between spin-up and spin-down electrons and thus creates 
an imbalance. From Maki's calculation(Eq.(23) in ref([T^|), we get that the 
coefficient of the fourth order term(in ip) the free energy expansion is given by 
— [(mp )/(87r 2 (27rT) 2 )]_Re^^L (n + \ + ip)~ 3 where m,po are the mass and 
Fermi momentum of the electron, respectively, and p = /iH/2ttT. Near the 
tricritical point, we can approximate the above form to — 0.002?7jpo after some 
algebra. In our theory, an expansion of the free energy in powers of ip 2 gives 
the coefficient of the fourth order term to be -(B 2 /2A 3 )(bA 3 /2B 2 - h 2 ) which 
can be approximated as — (bB 2 /2A 3 ) 2 near the tricritical point. So we can 

put the following condition on the phenomenological coefficients of our theory 
1 

(ll 3 ') = 0.002mpo- We want to mention that as Maki's calculation was based 
on superconducting system, so this estimation is valid only in the deep-BCS 
limit. More recently, Sheehy[li| had done a similar calculation. 

Finally we mention one important point regarding this free energy. As 
our proposed free energy does not contain any gradient terms, it will not re- 



produce the Fulde-Ferrell-Larkin-OvchinnikovfFFLO) [20i l21j superconducting 
state which is a superconducting state with a modulated order parameter. But 
one can extend this model free energy to include the gradient terms and study 
the FFLO state(22j, [23|, l24j. Presently we are working in this direction. 
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3. Conclusion 

To conclude, we have analysed the finite-temperature phase diagram of the 
two-component, imbalanced Fermi gas, as a function of population imbalance, 
using a phenomenological Ginzburg-Landau free energy. The novelty of our 
approach lies in using two order parameters, which, we assert, are actually 
embedded in the system: one superconducting and another imbalance order pa- 
rameter. This approach is inspired by the Blount and Varma's Ginzburg-Landau 
free energy which was proposed to describe superconducting-ferromagnetic tran- 
sition. Our analysis reproduces the phase diagram reported by experiment. We 
also point out that the specific heat behaves anomalously and provide a de- 
tailed analysis of the specific heat. We anticipate that the proposed free energy, 
with two competing order parameters, will lead to improvements in our current 
understanding about the imbalanced Fermi gases. 
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Appendix A. 

If we integrate out magnetization from our model free energy, we get an 
effective free energy which is different from the standard 0(ip 6 ) theory: 



e 



p me -(f<A 2 +!</> 4 +4m 2 -™'*+#™V) 



= const xe V 2 " 4 " J . (A.l) 

So the effective free energy reads 

„ a , o b , A h 2 , . s 

F 'ff=2^ + - 4 f- 2(A + B^Y (A - 2) 

The interesting point is that we can recover the standard 0(ip 6 ) free energy from 
this effective free energy by expanding F e f / near transition, where the value of 
ip is small: 

h 2 fa Bh 2 \ , 2 (b B 2 h 2 \ . B 3 h 2 iP 6 „ r(l8 , AO , 
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Now minimisation of the effective free energy with respect to ip can be done as 
follows: 



=*> ip[(a + bip 2 )(A + Bil} 2 ) 2 + Bh 2 ] = 0. (A.4) 

So, we have shown that the non-zero solutions of superfluid order parameter 
will be given by the solutions of eq. ([2]). 
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